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By assuming the existence of a pseudopotential smooth enough to do Hartree-Fock variations and 
good enough to describe nuclear structure, we construct mass formulae that rely on general scaling 
arguments and on a schematic reading of shell model calculations. Fits to 1751 known binding 
energies for N,Z> 8 lead to rms errors of 375 keV with 28 parameters. Tests of the extrapolation 
properties are passed successfully. The Bethe-Weizsacker formula is shown to be the asymptotic 
limit of the present one(s). The surface energy of nuclear matter turns out to be probably smaller 
than currently accepted. 
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The local structure of the mass surfaces is fairly smooth and amenable to algebraic analyses of shell model origin 
HU that lead to precise mass formulae - the rms errors may go below 400 keV (J] - but they need many parameters and 
have questionable extrapolation properties far from stability. The droplet (FRDM) || and Thomas-Fermi (EFTSI) 
H mass tables are designed to reach the neutron drip-line, as demanded by calculations of r-process nucleosynthesis, 
but they are not very precise: with 30 parameters the former yields an rms error of 673 keV for N, Z > 8 while the 
OO ! latter needs only 9 parameters for an rms error of 730 keV for A > 36 . 

In a recent fit it was shown that with 12 parameters the rms error for N > 28, Z > 20 could go down to 385 
keV. This fit does not extrapolate well and assumes a-priori knowledge of the right closures and the right boundaries 
between spherical and deformed regions. In a companion paper || rigorous microscopic guidelines were proposed to 
take advantage of the simplicity detected in [^|, while eliminating the drawbacks. A brief review of these guidelines 
, will make clear which is the problem to be solved to obtain a good mass formula. 

■ The starting point is a basic separation property of the interaction (anticipated in Q but proven in Q ): A 
(renormalizcd) Hamiltonian, ready for use in a shell model calculation can be written as H = H m + Hm ■ 

■ The monopole part, 

ON 



H m = V ]a k im k (mi - S k i) + b k i(T k .Ti - ^^-8 k i) (1) 



Sm k 



k,l 



contains only quadratic forms in the number (m k ) and isospin (T k ) operators for orbits k. The expectation value 
of H m for any state is the average energy of the configuration to which it belongs (a configuration is a set of states 
with fixed m and T for each orbit). In particular H m reproduces the exact energy of closed shells and single particle 
(or hole) states built on them, since for this set (which we call cs ± 1) each configuration contains a single member. 
We have emphasized the fact that eq.(l) is ready for shell model calculations because its form is not the usual one, 
in which effective core and single particle terms are separated from the two body parts acting in the valence space. 
It should be noted that the kinetic energy can always be written as two-body by eliminating the center of mass 
contribution. 

The multipolc Hamiltonian Hm contains all other terms (pairing, quadrupole,etc), and does not affect the cs ± 1 
states that bound the shell model spaces. 

This general result depends only on the assumption that there exists a pseudo-potential smooth enough to do 
Hartree-Fock (HF) variation. The renormalization process that defines the effective interaction leads also to many 
body forces that vanish at the closed shells. Their effects cannot be distinguished from those of the other higher rank 
terms that will appear later in simulating configuration mixing [||. 

In Q it was shown that Hm is well given by the realistic nucleon-nucleon interactions, and therefore parameter-free 
while H m must be treated phenomenologically because of the poor saturation properties of these forces. The problem 
is that eq.(l) contains hundreds of parameters, each of them a function of A and T (T — \N — Z\/2). In the context 
of U local fits are possible because only a fraction of the parameters are active and information on the cs ± 1 states 
is available. To say something general about H m and, in particular, to construct a mass table, we must proceed 
differently. 

In section I we shall find how to extract from eq.(l) the few dominant terms responsible for bulk properties and shell 
formation. In section II we show how to mock configuration mixing through simple algebraic expressions, and how to 
determine the spherical-deformed boundaries. Section III contains the results of a fit, whose asymptotic behaviour is 
discussed in section IV. 
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I The form of H m . 

To discover the form of the dominant terms in H m we rely on two geometric properties of the realistic forces. 

1. The A, T dependence. The aki and bki coefficients in eq.(l) must behave as typical matrix elements 

V(uj)klmn = —V(uj ) k i mn + 0(uj 2 ), (2) 

a result from ref. ||,but adding an O(oj 2 ) correction warranted for large oscillator constant u>. Then we know that 

^ = ^!.40A-V3 (3) 
< r 2 > 

a classical result Jll]] obtained with the standard < r 2 >= 0.86A 2 / 3 . A much better fit to the radii is obtained with 
< r 2 >= 0.90R C 2 where R c = A 1 / 3 (l - (T/A) 2 ) § . 

By combining this result with eqs.(2,3), and calling a typical term in eq.(l) TT(A,T), we obtain 

T(A,T) = (T/R)(l-p(T)/R)) , R = R^/A 1 ' 3 (4) 

for the amplitude T(A,T) affecting operator T. (The hat "will be omitted when no confusion is possible). 

2. Diagonal form. The isoscalar master terms. The matrices aki and bki in eq.(l) can be reduced to diagonal form. 
For ciki, say, we have 

^a k im k mi = ^ e Al (^m fc / fe)j ) 2 (5) 
k,l A 1 k 

and we borrow from |l(]] the result that the largest contribution (by far) to a realistic force has the form: 

t T (A) = (J2 m p /^/D;) 2 eo, (6) 
p 

where m p is the number operator for the oscillator shell of principal quantum number p ,and D p the degeneracy 

{D p = (p + l)(p + 2)). Calling M = m P /\/Dp~, setting nip — Hp + Zp, where n, z are number operators for 
neutrons (v) and protons (tt), and filling shells p v and p v up to some pj v and pf„ Fermi level, we find : 

< M >^ i[(3iV) 2/3 + (3Z) 2 / 3 ] = 0{A 2 ' Z ) (7) 



where we have approximated yj (p + l)(p + 2) ~ p + 3/2, and used N = Yl n p v — J2(p» + ^)(pv + 2) = |(p/„ + 2) 3 
and Z = \(p u +2) 3 . 

Therefore < f > = 0(A i / z ), and, from eq.(4) we have < f > T(A,T) = 0(A). As we shall see, this term is 
responsible for most of the bulk energy. At the same time it produces strong magicity at the harmonic oscillator 
(HO) closures through the \JD V scaling of the m p operators. This effect is absent in an infinite medium, where 
a ki = 0(A~ l ) and the total number operator m = replaces the eigenvector M to ensure the O(A) behaviour 

of the energy. 

If M is viewed as a coherent state, by symmetry the other eigenvectors that can be constructed out of m p operators 
are simply these same operators orthogonalized to M. Then we expect another monopole term of the form P = 
£ f{p)m p 2 /D p , which goes either as ^ 2/3 if f(p) = 0(1), or as A if f(p) = 0(p). 

By now we can state the general rule: 

Monopole terms must be symmetric quadratic forms of properly scaled operators. 

Proper scaling is determined by the assumed form of the amplitudes (i.e. eq.(4)), and by the asymptotic properties 
expected (or postulated) for each term. It is worth pointing out that proper scaling is a delicate matter that must 
be decided on empirical grounds, and to illustrate the possible ambiguities we consider a particular a k i coefficient in 
eq.(l) and take k and I to belong to the same major shell. As contributor to the main term M 2 (i.e. eq.(6)) a k i must 
contain a part that goes as u>/D p — 0(A^ 1 ) 1 but it also contains a contribution to P — ^ f(p)m p 2 /D p , that goes as 
f(p)Lu/D p , i.e. 0{A- r ) if f(p) = 0(1), or 0{A- 2 ^ 3 ) if f(p) = 0(p). (In ref. ||, an A~ - 75 scaling is extracted from 
experimental spectra, quite consistent with f(p) = 0(p)). 

Spin-Orbit and isovector terms. In deciding which are the contributions that are necessarily present beyond the 
ones already discovered, we consider eigenvectors orthogonal to M . To obtain the observed extruder- intruder closures 
(EI) we need spin-orbit (SO) effects and we introduce 
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Sp = pm.jp — 2m rp = (I ■ s) p . (8) 

Here (refer to Fig.l) jp is the largest orbit in the p-th shell and rp regroups all the others. (I ■ s) p is the operator 
that produces the same splittings as (/ • s) p and then collapses the r-orbits to their centroid value. The rationale for 
considering only two types of orbits is clear from Fig.l : we want to give top priority to shell formation (i.e. the cs 
part of the cs ± 1 set ).The combinations of operators other than nij and m r will contribute to subshcll effects 
that we incorporate in Hm and treat later. It should be noted the formation of the EI closures is not inconsistent 
with the survival of HO ones in the ligth nuclei. 

The other necessary ingredients are the isovector counterparts of m p and s p . Calling t p — \n p — z p \, 

S p = p{7~ijp -\- zjp) l 2{n r p -\- z r p) , and 

Stp p | tXjp Zjp | 2 1 TL rp Zrp I • 

we introduce the variables 



MT n 



ST n 



Stp 
2(p + l) 



(9b) 



Finally we construct the 12 possible symmetric quadratics shown in the first part of table I. The abreviations are 
(F=full,P=partial), (M=master, 5=SO, C=cross), (y4=isoscalar, T=isovector). The C terms are hybrids that could 
not represent eigensolutions in eq.(5), but they could modify (i.e. mix) slightly the M and S terms. They have been 
introduced to ensure - through FCA and FCT - the presence of representatives of the conventional spin-orbit force, 
including possible isospin effects. It is only for "symmetry " reasons in table I that PC A and PCT are also present, 
but their effects are not expected to be significant and they will excluded from the fits considered in this paper. 

FSA and FCA are scaled so to go as the ordinary (l-s term with ( = 0{A~ 2 / 3 ) Q. The FXT operators are 
scaled as the FX A ones. For the PXY contributions we introduce a factor to allow for possible ambiguities in 
the f{p) factor discussed above. The fits will decide in favour of a = 1/2. 

In addition to these terms, H m includes a pairing (V^,) and Coulomb (V c ) contributions as well as a 4T(T + 1) term 
whose presence we explain below. Although V p belongs in principle with Hm, it has been incorporated here because 
it affects equally all nuclei. Similarly V c contains multipole parts, but at the level of precission we work they are 
negligible. 

Asymtotic form of the master terms. Except for a small contribution from PMA and PMT (when a = 1/2), the 
asymptotic expressions 



(§^) 4/3 (l 



(§^) 4/3 (iS: 



§(2W) 



and 



FMA 

FMT : 

determine the bulk and symmetry energies of nuclear matter. As these quantities must emerge from a balance of 
potential and kinetic energies that our arguments 



TABLE I. The operators F (called T here) in H m , H s and H d 



FMA = (J2 MA v) 2 
PMA = J2(MA P ) 2 D™ 

FMT = (J^MTpf 

PMT = J2( mt p) 2d ? 
4T(T+ 1)A~ 2/3 



H s (n — D v — n, z — — z) 
SQ+ = 2(nn) 2 Df " 3 + 2{zz) 2 Df ~ 3 



H d (ri =n- JU, n' = n+ JU), 
QQl = {n'n'D- 3/2 ± z 'z ' D~ z/2 ) 2 



a = 1/2 

FSA = (J2SA P ) 2 
PSA = YiSArfD" 

FST = C£ st p) 2 
PST = J2(ST P ) 2 D« 

V p = -mod(N, 2) - mod(Z, 2) 

S2 = nnD' 1 + zzD^ 1 
SQ- = Annzz(D v D 1T ) l3 - 3/2 

(z' = z- JU, z' = z + JU) 

QQ± = QQ0±QQ1 



R c = A 1/3 (l - (T/A) 2 ) 

FCA = J2 MA p Dp 1/2 J2 SA p Dp 1/2 

PC A = Y J {MA P ){SA P )D^- 1 

FCT = MT P D~ 1/2 ST P D P 1/2 

PCT = X;(MT P )(ST P )D^ 1 

V c = [-Z(Z - 1) + .76[Z(Z - l)] 2/3 ]/i? c 

S3 = nn(n - n)Dl {p - 1] + zz(z - z)Dl {P ~ 1] 
= 1/2 

D3 = n'n'(n' - fi)DZ 2 + z'z'(z' ~ z')D- 2 
DO = 16, JU = 4. 
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have bypassed, the shell effects produced by FMA and FMT may be partially spurious and it is convenient to 
decouple them from the asymptotic contributions by adding these as independent terms. 

The three operators that turn out to be favoured by the fits are 4T(T + 1)/A 2 / 3 , and FM+ and PM+ where 
XM+ = XMA + XMT, {X = F, P), whose asymptotic form is 

XM+ = 1.717A 4 / 3 + .3816(2T) 2 /A 2/3 + 0(A 2/3 ) (10) 

II Configuration mixing due to Hm 

Once H m has ensured shell formation, mostly of EI type, the variables that become important in modelling config- 
uration mixing are n v and z v , the number of valence particles in EI spaces of degeneracy D v and (see Fig.l). For 
the light nuclei these model spaces are unconventional but not necessarily wrong: 12 C can be a good core [Q, and 
the N, Z=14 closures are at least as convincing as the N, Z=20 ones ||]. 

Spherical nuclei. H s and subshell effects. The energies of the spherical nuclei will be mocked by H s , a linear combina- 
tion of the four operators S2, S3, SQ+, SQ-, listed in the second part of table I, each affected by a coefficient as given 
in eq.(4). As explained in ||, H s is largely devoted to smooth subshell effects. Since S2 is easily absorbed in H m , 
and SQ+ is not expected (or found) to be important, we are left with only two terms whose main task is to correct 
the artificial filling patterns due to the imposed degeneracy of the r-orbits (we shall return to this point in section III 
and in the conclusion). The D 13 factors in table I reflect scaling uncertainties associated to these corrections. 

Deformed nuclei. H d . The energies of permanently deformed ground states will be described by H d whose form was 
not derived in || , but it can be easily shown that it must must be similar to that of H s , although the underlying physics 
is different. The argument is that the onset of rotational motion is associated with the interruption of normal spherical 
filling by the promotion of JU^ protons and JU V neutrons to configurations of the next HO shell. Nilsson diagrams 
indicate that JU^ =JU V =JU=4. The intruders bring a loss of monopole and a gain of quadupole energy represented 
by the constant DO, (the number 16 in table I is an arbitrary factor). The deformation energy of the particles in the 
lower orbits is simulated by quadrupole terms QQ + and QQ-, corresponding to the filling of equidistant Nilsson 
orbits, and D3 is in charge of the balance of monopole effects. Scaling is such that deformation energies have the 
standard A 1 / 3 behaviour. 

The fitting procedure and the spherical-deformed boundaries. Energies taken to be positive are given by the 
expectation values 

E(N, Z) =< H m > + <H s >(l- S d )+ <H d >S d 

= max(< H m > + < H s >,< H m > + < H d >) (11) 
TABLE II. Parameters of the 28p and 28p* fits / A => FM+ (28p) or A 4/3 (28p*) / V p and V c given in Table 3. The 



numerical values correspond to binding energies in MeV. 
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A 


PM+ 


4T(T+1) 


FS+ 


FC+ 


S3 


DO 


QQ- 


PS+ 


PS- 


FS- 


D3 


SQ- 


QQ+ 


T28 


9.55 


-0.77 


-37.23 


6.03 


-11.18 


0.47 


-38.1 


25.5 


-0.9 


-0.13 


1.4 


-0.9 


0.35 


4.6 


p(r)28 


0.89 


1.3 


1.38 


4.55 


5.11 


4.75 


4.81 


4.09 


5.24 


5.03 


4.21 





4.47 





T28* 


16.73 


-0.78 


-33.35 


6.05 


-18.04 


0.42 


-39.8 


13.6 


-1.2 


-0.16 


1.7 


0.3 


0.16 


3.3 


P (r)28* 


1.44 


4.87 


1.45 


5.39 


4.13 


4.34 


4.79 


2.75 


6.09 


5.30 


4.24 





4.35 
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the lowest possible orbits are filled for spherical nuclei (Sd = 0), while for deformed ones (Sd = 1), JU particles 
are promoted to orbits j. The calculations are conducted by initializing Sd, fitting E(N, Z) in the first equality of 
eq.(ll) to the 1751 mass values for N, Z > 8 in the latest compilation Q], then resetting Sd through the second part of 
eq.(ll) and iterating until convergence. Deformed states are only allowed in the regions where neutrons and protons 
fill different orbits i.e. above A w 100. Below, the onset of rotational motion is not sharp, pairing effects are stronger 
and a specific treatment would be called for. 

III The fit. It has already been mentioned that only a particular combination of the master terms was favoured, 
and that the 52, SQ+, PC A and PCT contributions would be omitted. In addition, we keep only one XCY term 
and omit the p factors in D3 and QQ+ (the less significant of all the included operators). This leaves us with 28 
parameters, shown in table II. This fit, 28p, yields an rms error of 375 keV. In the 28p* fit, FM+ is replaced by A 4 / 3 , 
leading to an rms error of 485 keV. The interest of comparing 28p and 28p* will become clear in section IV. 

The notation for the operators is XY± = XYA ± XYT and we have preset JU = 4, a = f3 = 1/2. 

It is clear that the p ratios fall in two broad categories: for the master terms they are close to unity, and for all 
others they bunch in the interval 4.5±1. According to eq.(4), it means that the corresponding contributions change 
sign at p ss A 1 / 3 , i.e. A w 100, the region where the j orbits can start filling before the r ones are full. This is a clear 
indication that parametrization (4) acts in a way unrelated to the derivation in section I, and that all contributions - 
in addition to doing what they were designed for - are busy simulating subshell effects. Setting p—0 is in general of 
little consequence , except for 53, which entails a loss of 300 keV, and to a lesser extent QQ— (loss of 150 keV); and 
it is possible to obtain fits with rms errors of some 600 keV, with only a dozen parameters. 

Some elements of evaluation are given next. 

Radii. The radius extracted from V c , r « 1.235i? c is close to the fitted r w 1.225i? c 0. The use of R c rather than 
A 1 / 3 leads to a gain of 34 keV. 

Deformed Nuclei. The number of deformed nuclei is 396 for 28p and 393 for 28p*, with rms errors of 254 and 313 
keV respectively. The sharp transitions (e.g. at N = 89,90) are accurately reproduced. 

Extrapolation properties When fitted to the 1503 nuclei in the 1983 compilation, 28p yields the rms error of 357 
keV. When the same parameter are used on the 248 new nuclei in the 1993 tables the rms error is 470 keV, which 
amounts to good predictive power by present standards. The increase in error is due only to spherical nuclei, and 
we have found that discrepancies exceeding 700 keV are almost invariably associated with subshell effects that are 
important only in spherical regions. Progress in this matter can be expected through a more sophisticated treatment 
of the 5 operators. 

When it comes to systematic behaviour, which includes many spherical nuclei and all the deformed ones (unsensitive 
to subshell details), the fits do very well as emphasized by another test: when 28p in table II is used on the 2542 
entries (i.e. including systematic trends) of the 1993 tables the rms error is almost unchanged (407 keV). 

IV Asymptotics and the Liquid-drop. 

The 28p and 28p* entries in Table III are the asymptotic contributions, given by eq.(10) for FM+ and PM+, plus 
the Coulomb and pairing terms. Under 6p we have the results of a direct fit to the data by varying these 6 



TABLE III. Asymptotic forms of the fits compared with a pure LD form (6p). (T2 is for 4T(T+1)) 





A i/3 /R 


A 4/3 /R 2 


T2/A 2/3 R 


T2/A 2/3 R 2 


V p 


v c 


28p 


15.07 


-12.79 


-33.88 


48.51 


5.18 


.699 


28p* 


15.39 


-17.54 


-33.64 


49.81 


5.15 


.699 


6p 


15.42 


-17.56 


-33.65 


50.91 


5.18 


.699 
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parameters, that can be identified to the classical liquid drop (LD) ones by expanding the R and R 2 denominators 
(rms error of the 6p ht: 2499 keV). 

The agreement between 28p, 28p* and 6p in Table 111 is excellent - even spectacular - except for the surface 
coefficient in 28p. (Replacing FM+ by FMA leads to an intermediate situation with an rms error of 436 keV and a 
surface coefficient of -14.6). 

It is apparent that conventional LD parameters as extracted in 6p are fully reproduced by 28p*, a very good fit 
by present standards. However, 28p is significantly better and provides a strong hint that the surface energy may be 
smaller than the conventional value. 

Conclusions and prospects. 

The rms errors we have obtained may seem good - even impressive - by present standards, but they are still large 
compared to those obtained in shell model calculations, i.e. some 300 keV in the lightest nuclei and 150 keV around 
A = 100 0. The reason is that our treatment of H m is too crude since we have forced an artificial degeneracy on the 
r-orbits. As explained in ref. the correct behaviour can be simulated to some extent by the H s and Hd terms - 
designed to describe configuration mixing - but also capable of coping with "monopole drift" , i.e. the smooth part of 
subshell effects. 

A more refined treatment of the S p operators in H m will eliminate the large p ratios, which are physical - in the 
sense that they represent a real effect- but mask the true contribution of a given term, and introduce unnecessary 
phenomenology in the formulation. 

It should be clear from these remarks that the fit(s) we have proposed are only exploratory. There are grounds to 
expect much better ones. 

NOTE. It takes only seconds to generate a mass table. The program is available from dufloj@frcpnll.in2p3.fr . 
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HO EI 

r(p+l) 

D p+1 ^L±1^ 

j(P+l) 



r(p) D v ^ 
j(p) 



D jp = 2(p+l) D rp =p(p+l) 

FIG. 1. HO and EI major shells. 
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